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Abstract 

Euler integrals of deterministic functions have recently been shown to have a wide variety of 
possible applications, including in signal processing, data aggregation and network sensing. 
Adding random noise to these scenarios, as is natural in the majority of applications, leads 
to a need for statistical analysis, the first step of which requires asymptotic distribution 
results for estimators. The first such result is provided in this paper, as a central limit 
theorem for the Euler integral of pure, Gaussian, noise fields. 
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1. Introduction 

The Euler characteristic x(^) of ^ nice set A is perhaps the oldest, and most funda¬ 
mental, of its topological invariants. For a compact A C M^, the Euler characteristic is 
merely the number of its connected components (each one of which will be an interval, pos¬ 
sibly containing only a single point). For A C M^, x(^) becomes the number of connected 
components minus the number of holes, while in three dimensions x(^) can be written as 
the alternating sum of the numbers of components, handles and hollows. Similar (and, of 
course, more precise) definitions as alternating sums of Betti numbers, numbers of facets 
of simplices of differing dimension (when A is triangulisable) or as indices of critical points 
when a Morse theoretic setting is appropriate, extend the Euler characteristic to a wide 
variety of sets in arbitrary dimensions. 

However, more important for us is that the Euler characteristic is also a valuation, 
which means that, when all terms are defined, we have the additivity property, 

x{AV3B) = x{A) + x{B) (1.1) 
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Given additivity, it is natural to attempt to use x to define an integral on a suitable family 
of functions, and, indeed, to a large extent this can be done. The resulting theory is known 
as Euler integration. 


1.1. Euler integration 

Although in many ways Euler integration has its roots in classical Integral Geometry, 
a more complete and modern theory began to evolve in the 1970’s. More importantly for 
us, however, is that it has experienced a rapid development in the past decade from both 
applied and theoretical aspects, providing for some ele gant and novel results. We shal l not 
attempt to survey these here, since the recent papers of Baryshnikov and Ghrist ( 2009l l and 
Gurry et al.l (120121 ) provide excellent and broad expositions. Rather, we shall go directly 


to two definitions. 


Definition 1.1. Let M C M” he compact, with finite Euler characteristic. Then a con¬ 
tinuous function f : M — )• M is called tame if the homotopy types of f~^{{—oo,u]) and 
f~^{[u,oo)) change only finitely many times as u varies o?;er M, and the Euler character¬ 
istic of each set is finite. 

Definition 1.2. // / : M — )• M is tame, then the upper Euler integral of f over M is 
defined by 

[ f\dx] = f [xif >u)-xif <-u)]du, ( 1 . 2 ) 

J M JU=0 

where 


Xif<u) = xif ((-oo,u])) 


and 


Xif>u) = xiM)-x{f<u). 

Reading in between the lines that do not appear in the above definition, one would guess 
that there is also a lower Euler integral (there is!) and that there has to be a more direct way 
to define an integral that follows from the additivity of (|l.ll) . In fact, this is also true, and, 
as a result, the Euler integral shares many common properties with the classical theories of 
integration. However, it is somewhat more delicate, since although (jl.ll) extends to a finite 
inclusion-exclusion form, it does not typically extend to the countably infinite case needed 
for a standard measure based theory of integration. The definiti on that we have chosen 
above avoids these issues, and in taking it we follow the lead of iBobrowski and Borman 
(|2012l l who, by taking (II.2p as a definition rather than a property, save often irritating but 
unimportant (for our needs) technicalities. 
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1.2. A motivating application 


An inte resting appli c ation of the Euler integral is described in iBarvshnikov and Christ 
( 2009! ) and Curry et al. ( 2012 ). 

Suppose that an unknown number of targets are located in a region M C M"", and each 
target a is represented by its support Ua C M. Suppose also that the space M is covered 
with sensors, reporting only the number of targets each one sees. Let /i : A —)• Z be the 
sensor field., i.e. 


h{x) = {targets activating the sensor located at x} 


Then, if all the target supports satisfy x{Ua) = for some ,8 7^ 0, the readings from all 
the sensors can be combined to obtain the exact number of targets via the relationship 


N = {targets] = ^ [ h\dx\. 

P Jm 


(1.3) 


Note that we do not need to assume anything about the targets other than that they all 
have the same Euler characteristic. Eor example, we need not assume that they are all 

convex or even have the same number of conne cted components. _ 

While everything in (11.31) is deterministic, Bobrowski and BormanI (2012) raises the 
question as to what happens when the deterministic ‘signal’ x = h\dx ], is observed via 
a noisy measurement Y = jj^{h + X)\dx], where A is a smooth random process on M. 
They show that, although Euler integrals are not always additive, in this case it is true 
that 

Y = [ {h + X)\dx] = [ h\dx]+ [ f\dx] = /3A + C, 

Jm Jm Jm 

which leads to the obvious estimator A of A given by 


A = /3 


-1 


y -E 


Um 


x\dx\ 


(1.4) 


The main result of Bobrowski and Borman ( 2012l i is an elegant calculation of the ex¬ 
pectation in (II.4p when A is a smooth Gaussian or Gaussian rela ted random held and M a 
stratihed manifold, based on the Gaussian kinematic formula of I Adler and Tavlorl (j2007l L 
Their computation leads to an explicit, closed form (and often somewhat surprising) ex¬ 
pression for the expectation. We shall have more to say about this in Section 01 

Motivated by the above, what this paper concentrates on is a central limit theorem 
(henceforth CLT) needed to go from the estimation provided by A to inference. 
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1.3. A CLT for the Euler integral 

The main result of the paper is formulated in Theorem 13.11 which states that if X is a 
real valued, almost surely stationary Gaussian random field on M”, satisfying certain 
technical regularity and decay of memory conditions, and if we define 

T[o,h„[X] ^ [ 

J [0,m]' 

then, as m —oo, 

rnn/2 

for some limiting variance a\ > 0, where ^ denotes convergence in distribution. 


X\dx\, 

^ iV(0,4), (1.5) 


l.f. Outline of the paper 

Before giving the formal version of (jl.5p in Section [Sj in the following section we shall 
set up considerable preliminary material, treating the regularity conditions that we require 
on X as well as background material on Wiener chaos expansions and the Morse theo- 


Houdre and Perez-Abreu ( 

I994I') and 

Mai or 

(2014j) 


to develop the chaos expansion for the 
upper Euler integral of a Gaussian random field. This, to gether with some regularity and 
convergence results, are combined with a general CLT of Nourdin and Peccati ( 20121 ) for 
chaos expansion s, to make up the proof our main CLT. Many o f the proofs here owe a lot 
to the p apers by Kratz and Leo J ( IOOTI ). Kratz and LeonI ( 2001 ). and especially the recent 
work of Estrade and Leon ( 20151 ) . 

In Section 0] we look at a direct calculation of the mean value of Euler integral (for the 
isotropic case), showing its depende nce on the order-one Lipshitz-K illing curvature of M, 
and discuss the surprising results of Bobrowski and Borman ( 2012l f alluded to above. 


1.5. Acknowledgements 

1. A reader familiar with the paper Estrade and Leon ( 2015l l will note a strong re¬ 
semblance between many of the technical parts of that paper and those in Section [3] in 
the present paper. This is not accidental, since that paper develops a CLT for the Euler 
charaeteristic of excursion sets, while we treat a CLT for the Euler integral. Neither result 
implies the other, and, as far as we can tell, neither is derivable from the other without a 
lot of additional work. Nev ertheless, many of t he de tails of the calculations are similar. 

We became aware that Estrade and Leon ( 20151 ) was in preparation early in our own 
work, and are very grateful to Anne Estrade and Jose Leon not only for sending us a 
preprint of the final paper, but also for sharing in-progress versions long before the final 
preprint was ready. Doing so not only saved duplication of effort, but made our lives much 
easier. 

2. Research supported in part by FP7-ICT-318493-STREP and ERC 2012 Advanced 
Grant 20120216. 


4 































































2. Gaussian random fields and chaos expansions 


Before we set up our results in a formal fashion, we need some preliminaries. In partic¬ 
ular, we require a collection of regularity conditions on our random fields which will make 
the Euler integral well-defined, amenable to analysis, and which are sufficient for our CLT 
to hold. 

In addition, and this will take up most of the section, we need to set up a number 
of results related to chaos expansions. These will be used in the remainder of the paper 
to express the Euler i ntegr al in this form and then prove our CLT via a general CLT of 


Nourdin and Peccatil (j2012l ') for chaos expansions. 


For general preliminaries on random fields and their co n nectio n to Morse theory we shall 
use th e often complementary books bv lAdler and Tavlorl (j2007l l and lAzai's and Wschebor 
( 20091 ). while for a good treatment of the Wiener chaos we rely on Nualart ( 200f)l ). Besults 
below that we refer to as “standard”, “well known”, or for which we fail to offer even these 
descriptions, can be found in one of these references. 


2.1. Tame Gaussian fields 

For the remainder of this paper, X will denote a real valued, mean zero, unit variance, 
Gaussian random field on M”’, n > 1. We denote its covariance (and correlation) function 
by p : MT —M. For a function / : M” —)• M we denote its gradient by V/, writing this and 
other vectors as row vectors, and its Hessian by V^/. We shall occasionally treat V^/ as 
a vector rather than a matrix, in which case, because of symmetry, it will have n{n + l)/2 
elements. It should be clear from the context whether we are using the matrix or vector 
interpretations. Generic constants, which may change from line to line, are denoted by C. 

We write Gov(y) for the covariance matrix of a random vector T, and the ubiquitous 
symbol | • | to denote all of modulus (of a real number), length (of a vector) and determinant 
(of a matrix). Again, usage should be clear from the context. 

The regularity conditions we shall require on X are summarised in the following defi¬ 
nition. 

Definition 2.1. Let X = {Ai(t), t G M"'} he as above. Then we call X tame if the 
following conditions all hold. 


(i) At each t G the joint distribution of the vector {X{t), XX{t), V^A(t)) is non¬ 
degenerate. 

(a) The covariance function, p, of X is four times differentiable, and for some a > 0, 
and t small enough, each of its four-order derivatives satisfies 


p(4)(o)_pW(t) 


(4)/ 


< 


c 


(- In |t|)^+"' 


( 2 . 1 ) 
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(in) Set 




A 


sup 

0<m<4 


d^p 

■ ■ ■ ti^ 



( 2 . 2 ) 


Then ijj € and ipit) 0 as |t| ^ oo. (iv) 

(iv) Let NyCVX, M) be the number of points, t G M, for which XX(t) = v. Then, for 
any v G M"", 


E 


iN,iXX,M)Y 


< oo. 


There are a number of immediate, standard, consequences to tameness for a Gaussian 
random field. In particular, (ii) ensures that the trajectories of X are almost surely (hence¬ 
forth a.s.) in C'^(M”'), and, via the exponential integrability of the suprema of Gaussian 
processes (assured by the Borel-Tsirelson-Ibragimov-Sudakov inequality) that 

E[|supXi|^] < oo, 
teM 


for any compact domain M C MT, and any fc > 1. 

Condition (i) ensures that the realisations of X are a.s. Morse functions. We shall prove 
later that (in) ensures the decay of correlation necessary for a CLT to hold. Condition (iv) 
can be directly verified for specific covariance fu nctions using stan d ard in tegral expressions 
for the factorial momen ts of N„(XX, M ), as in lAdler and Tavlod ([2007|) [Theorem 11.5.1]. 
Alternatively, following Belvaevl (196^), this condition can be substituted by requiring 
non-degeneracy and smoothness for higher or der derivatives of the field . (See also a second 
moment calculation in an isotropic setting in Estrade and Leon (201^) [Proposition l.lj.) 


2.2. Correlation structure of stationary, tame, Gaussian fields 

In what follows, we shall often need details about the distribution of the random vector 


A ^ ((VA)i,...,(VA)„,A,(v2a)i,i,(v2a)i,2,...,(V2a)„,„)), 


of length Nn = 1 -|- n -|- n(n -|- l)/2, where 


(VA)i 


dX 

dti ’ 


and {X‘^X)i^j 


02 X 

dtitj ’ 


i,j = 1 ,... ,n. 


are the first and second derivatives of X. Since X is tame, all of these derivatives exist, 
and all are Gaussian. Furthermore, by stationarity, the distribution of Xg is independent 
of s, and the elements of the covariance matrix are given by derivatives of the covariance 
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function p at the origin. It is then well established that the covariance matrix A of X 
factorizes as 


A 



(2.3) 


where A(i) is the covariance matrix of VX and A( 2 ) is the covariance matrix of {X, V^X). 
Now let be a square root of A, and define the random field Y by 


y(s) = (2.4) 

The representation (12.31) induces a similar factorization on and Y. It is 

important to note that K is a vector valued random field. Furthermore, since, for each s, 
Y (s) is a vector of independent, standard normal variables, we shall call Y the decorrelated 
version of X. However, note that despite the independence of the elements of Y(s) for 
each s, the vectors H(s) and Y{t) are not independent for s 

For later needs, note that if we define the covariance matrix K by 

K{t) = Cov(y)(t) = E[yo'i"i], (2.5) 

then it is easy to check that its entries are bounded by 

< n2||A-(i/2)||2.^(i) 

< Cil){t). ( 2 . 6 ) 

Here ijj is given by (12.2h and (7 is a constant dependent only on the derivatives of p at the 
origin. 

Throughout this work, Y denotes the decorrelated version of X. Before concluding 
this section and while the definition of Y is still fresh in our memory, we introduce a new 
notation, and using this new notation state some of the relations between Y and X that 
we will need later on: 

• For an arbitrary vector u of dimension d and a set of indexes I = {ij}j^i where 
ij G {1,... ,n} and k < d, we define the vector Vx (u) by 

Vx{u) = (uii,tii 2 ,... ,UiJ. (2.7) 


In particular, with I = {N^ — n} we have 

V = U (Agf’rpi), 

and with I = {m, /} 


((VX)^,(VX)0 = Hx(A[;f¥(i)). 
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• For an arbitrary vector u and a set of indexes X = we implicitly define a 

symmetric matrix Mj (ft) constructed from the elements of u so that, when u = 
and I = {ij ,... ,4}, we have 


Mx (V(2)) = 


(V2A)i3,, (V^A) 


^2^2 




In particular, if X = {1,..., n}, then 


V^A = Mx(A[Jf)y(2) 


(V2A),,,A 


( 2 . 8 ) 


(2.9) 


^.5. T/ie spectral distribution ofY 

Since A, and so its decorrelated related version Y, are stationary, both have spectral 
representations. Whi le all pr o pertie s of the spectral representation of A follow from the 


classical theory (e.g. lYaglomI (jl962l ll we need to work a little to set up an appropriate 


representation for the vector valued field Y. In particular, we shall do this in the language 
of isonormal processes, which provides the necessary structure for later proofs. 

We start by noting that since A is tame, the function ip is integrable, and so by (12.6p 
the same is true of the covariance K of Y. Consequently, by standard spectral theory, Y 
has a matrix valued spectral density function / for which 


(A(r)),fc = E[YjiO)YkiT)] = [ e'<"-^>/,fc(A)dA, 


r G 


( 2 . 10 ) 


It is not too difficult to express the fij in terms of the spectral density of A (which is, 
essentially, the only ‘free parameter’ in the entire setup) but, fortunately, their explicit 
form will not be important in what follows. What is important, however, and follows for 
the non-degeneracy condition (i) of tameness, is that, for all A G M", {fjk{^)^k=i ^ 

symmetric, positive semi-definite matrix, and so has a symmetric square root (6jfc(-^))^fc=i- 
Consequently, we also have that 


= Y.^ji{X)hk{X). 

i 


( 2 . 11 ) 


2.4- Representing Y via the isonormal process 

Retaining the notation of the previous subsection, we start with a separable Hilbert 
space 4 of Hermitian functions 


h(j,A):{l,...,A„}xM-^C 


Hj, A) = h{j, -A), \\h\\sj < oo ^ , 


( 2 . 12 ) 









with the inner product 



(2.13) 


Next, we let j = I,... ,Nn a sequence of independent, real-valued, Gaussian white 

noises on and use them to define a random process over h £ Sj hy 



(2.14) 


the integrals here all being standard stochastic integrals. 

By construction W{h) is centered Gaussian and K[W{h)W{g)] = {f,g)s--j- Moreover, 
since the functions h are Hermitian the resulting random process is real valued. It is 
known as the isonormal Gaussian process on S). 

Finally, we want to relate Y to VF, as promised. To this end define a new family of 
functions, A), A; = 1,..., Nn, t G M”, in via 


<ft,m = A) = 6j^ra J, m = 1 ,... , iV„, A € M”. 


(2.15) 


It is straightforward to check that 




An immediate consequence of this is that the vector valued random field Y has the following 
particularly useful L 2 representation in terms of the isonormal process and the family (ps,k- 


Yi{s) ^ W{ips,l), l = l,...Nn, S£R^. 


(2.16) 


Note that it also follows from these calculations that 


= iE[yfc(o)yfc(o)] = 1, 


(2.17) 


and 




(2.18) 


where 6k^m is the Kronecker delta. 


9 



2.5. Operations on f € S) 

We now describe the basic operations on Sj that we will need later. Let {ej}j>i be an 

orthonormal family of functions in Sj, and write M = {1,... , Nn} X M”. Then, ej{X) : M —)■ 
C, j > 1, and we define the following operations: 

• Tensor product, ej (g) : M X M —^ C: 

[cj (g) efc](Ai, A 2 ) = ej(Ai)efc(A2), Ai,A2€M, (2.19) 

where ej (gi belongs to the Hilbert space with the inner product induced, 

component-wise, by the inner product in fj: 

(cj (g) e/c, e; (g) — {ej,ei')sj{efi,6jTi')sj ■ (2.20) 


• In a similar fashion, we define the m-fold tensor product of ej with itself: 

g.0...0g. ^ ( 2 . 21 ) 

m times 

where belongs to the Hilbert space with the inner product defined component¬ 
wise by the inner product in S^. Likewise, for G i^®^, G io®^, the tensor 
product of higher order in is 

(g) = e)- (g) • • • (g) ej (g) Cfc (g) • • • (g) efc . (2.22) 

q times P times 


• Take 0 < r < p < m. The r-contraction, (g) • • • (g) ej^) (8>r (eki <S> ■ ■ ■ <8> e^^), is in 

^0p+m-2r j, _ Jg ag 

(Cji (g) ■ ■ ■ (g) ejp) (g)o (4i ® ® ek^) = (Cji (g) ■ ■ ■ <8> ej^ <S> <S> ■ ■ ■ <S> 4^), (2.23) 


while, for 1 < r < p < m. 


(cji (g) • • • (g) ej ) (g)r (efc 


- ^ A 


, eki)sj 


U=i 


(g) • • • (g) ejp (g) 4,+i <gi • • • ® 4^). (2.24) 


For r = p = m we have 


(ej^ (g) • • • (g) ej ) (g)p (4i ® ® 4^) = 


Ski)s^ 


L«=i 


(2.25) 
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• Symmetrization of / G 

f = symm(/) = ^ ^/(A^(i), • • • , A^(q)), (2.26) 

where aq is all the permutations over the indexes {1, • • • ,q}. We write C 
for the space of all symmetric / € 

2.6. Wiener chaos expansion 

Take W{h) to be an isonormal Gaussian process on separable Hilbert space f). Write 

G = <j{W{h)), for the n-field generated by the random variables {W{h),h G ii} and 
for the space of all square integrable mappings from to M. 

We make use of {Hn}n£N^ the probabilistic Hermite polynomials, dehned by 

4.2 ^2^ 

exp(fx - —) = ^ —, X G M”, (2.27) 

2 m\ 

m=0 

and define 

OO 

Hci{x) = WHa^{xj), Xj G M”, flj G N, j = 1,... ,oo. (2.28) 

j=i 

Here, the sequence a = { 01 , 02 ,...} is such that only a finite number of elements differs 
from zero. One can make use of {Hg{x)} to construct an orthonormal basis for L^(^,M). 
The basis is given by the random variables 

OO OO 

<ha = H-,{W{ej)), o! = > (2-29) 

j=l i=0 

with W{h),h £ S) the isonormal process as dehned in Section ITTI and {ej}j>i the orthonor¬ 
mal basis of Sj. 

Moreover, the space L^(^,M) may be represented as the decomposition of a countable 
set of closed orthogonal subspaces {'Hm}meN, 

L\G,R) = e^^o'»-m (2.30) 

such that, for each m, {$a}|a|=m is a complete orthogonal system in Tfm C L^(^,M). Thus, 
given F G L‘^{G,R) there is a unique decomposition 

OO OO 

^ = E Z |a| = , (2.31) 

m=0 \ d\=m i=0 
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with ^ orthogonal for different values of m. 

We are now in a position to set up the Wiener chaos expansion. To this end define 
4 by 

<? 

Im{ej^ «) • • • «) CjJ = (2.32) 

i=l 

where ai = ^{A: : jk = i}- Using {/g(-)}, one can rewrite (12.3111 as 

OO 

f = (2.33) 

q=0 


This representation of the random functional F via the set of kernels {fq} through the 
family of linear operations, Iq, is called the Wiener chaos decomposition of F, and the 
operator Iq is called the multiple Wiener integral of order q. 

If the underlying Hilbert space is a Polish space of the form L'^{A,A,p), then Iq can 
be identified with multiple stochastic integrals. More specihcally, take {Ai} G A disjoint 
sets, and dehne 


Uq = symm 



.,iq 1 Al ® 


• ® lAq 


(2.34) 


The functions Uq are dense in symm[(L^(T, M,/x))®'^], and the integral is constructed hrst 
on the functions Uq by defining 


I{Uq) =Y,<^n,..A,W{A,)---W{Aq), 


(2.35) 


and then extending to a linear continuous operator on 


all of (L^(H,M,/x))®'^. We write 


I{fq)= [ ■ [ fq{tu--- ,tq)W{d^^)■■■W{d^^) U G T, (2.36) 

Ja Ja 


with W{dp) a Gaussian /x-noise. 

Below we list some o f the bas i c pro perties of Iq that will be of particular interest to 
us. For more details , see (INualartl. l2006l. Ch. 1) . For an in depth treatment related to the 
approach above, see Nourdin and Peccati ( 20121 ). 


• / G Imif) = Imif). 

• For any set of hi G such that ||hj||f 3 = l,Vi, we have 


HHaMihi)) = /|s|(symm(/xf“^ ® ® h®“-)). (2.37) 

i=l 
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• ^Ip{f)Iq{9) = 0 when p ^ q and EIp{f)Ig{g) = pl{f,'^sj'»p when p = q. 
Other properties of the multiple Wiener integral will be recalled when needed. 


It is in this language of Weiner chaos that we seek to represent the Euler integral, and 
then make use of the following result. 


Theorem 2.1 (Theorem 6.3.1. iNourdin and Peccatil ((2013)). Let {Fm)m>i be a sequence 
in L?‘{LL,Q,P) such that = 0 Vm. Suppose that the chaos expansion of F is given 

by Fm = suppose in addition that 


(a) -5- o-q as oo, for some > 0. 

(b) < oo- 


(c) 'dq>2 andr = l,---q-l, ||/™ /™|||®2(9-r) 0 as m oo. 

(d) lim sup^>i E^q +1 ^ 0- 


Then, 


Fm Af(0, Up) 


as m —)■ oo. 


2. 7. A Rice type formula and a Morse theoretical representation of the Euler integral 
Although up until now we have approached the Euler integral via integration theory, 
for the proofs to follow we require a slightly different approach, via stratified Morse theory. 
This theory links the topology of sets to the study of critical points defined on them. In 
particular, the Euler characteristic of excursion sets of the form 

A{M,u) = {f € M : f{t) > u}, 


is easily computed via properties of the critical points of /. 

To see how this works in our setting, we return to Section Fl. 11 with M = T„ = [0, m]"', 
and, noting that x(^n) = 1; obtain 




[1 - x{Tn n (/ 


< u)) - x{Tn n (/ < 


—u))]du. 


(2.38) 


We introduce p{s), the Morse index of a critical point s of /. (i.e. V/(s) = 0 and the 
Hessian V'^f{s) has p{s) negative eigenvalues.) We write T° for the interior of the cube, 
then proceed to decompose the boundary of T„ into open faces each of which is an open 
cube of dimension less then n. (The vertexes are considered to be zero dimensional closed 
cubes). To save on notation we denote any such face by J, and write {J} to denote the 
collection of all such faces. When quantities are evaluated with respect to a particular 
face, this will be denoted by an appropriate subscript. With the above notation, a careful 
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applic ation of Theorem 9.3.5 of Adler and Tavlor ( 200?! ) (see also Bobrowski and Borman 
(2012)) yields 


/ /rxi= E + 


{^■*eT°:V/(s)=0} 


+ E 1 E b (2.39) 

JG{J}\T° Y{sGJ:V/|j(^=0} j 

In the above, {tIjI are constant vectors at tached to every face of the cube Tn (the 
details can be found in Adler and Tavloil ( 200?! )). We identify 


{sGT°:V/(sl=0} 

as the contribution of the internal critical points to the Euler integral, and 

( 

E 

Je{J}\T° 


(2.40) 




(2.41) 


\{sGJ:V/|j(s)=0} 


as the contribution of the critical points on the boundary. 

From a critical point representation of this kind, one can develop an integral represen¬ 
tation of Rice type, and it is this that will be at the core of all the proofs to follow. 

Lemma 2.2. Let f be a Morse funetion. Then 

[_ f\dx] = lim + 




^6{J}\T° 


idin.jxdimj('^/|j(^)det(vV|j(s))/|j(^l{(v/(s),*)>o}rfs, 

(2.42) 


where is a {dim. J)-dimensional centered Gaussian kernel with the eovari- 

ance matrix u^Idim Jxdim j- 

Proof. Standard techniques for the construction of Rice type integral formulae, along with 
the fact that the cr —)• 0 limit of (pa'^i is a Dirac delta, establish that, for any J G {T}, 

lim/ (/><^ 2 i(V/(^) |det(VV(^)|/(^ds = / {6 oV f){s)f{s)ds 


(T—>-0 


E 

{5eT°:V/(s)=0} 


(2.43) 
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and 


lim / 0^2; 


(j —>^0 


•dim J X dim J 




— j ° ^/|j)(^/|j(^l{(V/(s),»7j)>0}'^^ 

= X] /(^l{(V/(s),rfj>>0}- 
{seJ;V/|j(s)=0} 


(2.44) 


Using the fact that the determinant of a matrix equals the product of its eigenvalues, it 
follows that signjdet (V^/|j) (^} = (—Thus, we can drop the absolute value in 
()2.43p and (I2.44p . and, applying p2.39p . complete the proof. □ 

We now have all that we need to formulate, and to prove, the main result of this paper. 


3. A CLT for the Euler integral 

Theorem 3.1. Let X = {APs)!.? € M”} he a tame Gaussian field, as in Definition \2.1\ 
Then, the (upper) Euler integral 

'h[o,™].[A] ^ [ X{s)\dx\ 

J [0,m]^ 

satisfies the central limit theorem 




m 


n/2 


W(0, (j|), as m —)• oo. 


where a\ > t) is defined by (|3.36l) below. 

Note that an expression for th e mean value of Euler integral, [A]], was derived 

Bobrowski and Borman ( 20121 ). and is also discussed in Section [4] below from the point 


of view of chaos expansions. 

Before starting the proof of Theorem 13.11 note that while expressions like (I2.39h and 
()2.42p relate to the full Euler integral, only the first sum in (I2.39P and the first integral in 
(I2.42p . which relate to contributions from the interior of T„, are relevant for the CLT. The 
reason for this lies in the normalisation of which applies equally to all terms. It 

follows from the calculations of this subsection that all non-interior terms, when normalised 
by converge in probability to zero, and so not affect the limiting distribution. We 

leave the (simple) details of this to the reader, and so in dealing with the CLT henceforth 
concentrate only on the interior terms. 

We start with a sequence of lemmas, which will ultimately be combined to provide a 
full proof of Theorem 13.11 in the following subsection. 
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3.1. Four supporting lemmas 

Lemma 3.1. Let X be a tame Gaussian random field. Let 


F(,,,mAX]= Y. {-^Y^^X{s), 

{s&T°-.VX{s)=d} 

and 

F(0m)4X] = [ <PMXX{s})detiV^Xis))Xis)ds. 

.1(0,m)" 

Then, 


(3.1) 


(3.2) 


F(o,m)^ W Y F(^Q jn)n [X] as fj —)• 0. (3.3) 

Proof. We deduce the L^ convergence from the following two facts: 

(^) F^o,m)YX\ —^ Yo„^^n[X] as (T —^ 0, 

(b) There exists e > 0 so that supg.E[F((j^^„[X]]^'''^ < oo. 


To prove (a), note that the trajectory of X is almost surely Morse and the result then 
follows from Lemma 12.21 To show (b), we write an upper bound for 
i ndep endent of a. To this end, using Federer’s coarea formula (cf. lAzai's and Wschebor 


( 2 OO 9 I I. Proposition 6.1, for a version couched in our terminology) we have 


' (0,m)" 


((>^2i(VX(s))|det(V2x(s))|ds = / 0^2i(n)iV^(VX,(O,m)’^)du. 


(3.4) 


Write Xsup = sup^g(o,m)" l^(s)|- Then, applying (13.41) . we obtain 

^ ((),2i(VX(s))det(V2x(sl)X(s)ds< |(Xsup)| / (/),2i(u)iV^(VX, (0, m)’^)du 

JR^ 


'(0,m)” 


Using the fact that a6 < ^ ^ for a, 6 > 0 and 4 + 4 = 1, we have 




sup; 


fMu)Nit{XX, {0,mr)du 

1 


(2+e) 


< - 
P 


(X,,p)(P(2+e)) 


H— 

q 


[ fMmu{XX,{0,mr)du 

JR" 


q(2+£) 


(3.5) 
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Taking expectations yields 


E 




(^^ 2 i(VX(s)) det(V^X(s))X(^ds 


(2+e) 


< + -A f (paniu)N^iVX, iO,mr)di 

P Q I JR" 


<?(2+£) 


(3.6) 


is finite due to our assumptions of tameness on X, so we focus on the 
second term in (13.61) . viz. 


E 


/ cPMu)N^i^X,{0,mr)di 

JR" 


q{2+6) 


(3.7) 


Jensen’s inequality, when applied to the inner integral (and not to the expectation), implies 
that the above can be bounded by 


E 


/ (/.,2i(u) [N^iVX, (0, mr)] . 

.JR" 

Using Tonelli’s theorem, this equals 

/ 0,2i(u)E[X„-(VX, (0,m)")]''^'+"^du, 

JR" 

and, finally (for q = 1 + e) under the assumptions of a tameness we have 

EX^;(VX, (0,m)”)(2+^) < M, 


and we are done. 


For the other faces J € {J}\r°, of dimension d = dim(J) < n, we have 


E 


<(’(t2i(VX(s)) det(V^X(s))X(s)l|i;vj(^(^) 


- (2+e) 


< E 




sup 


'/ 

JR‘ 


(j)^ 2 i{u)Nf^{VX, J)du 


(2+e) 


and can then repeat the same argument as above. 

For the next lemma, which deals with the Wiener chaos decomposition of T(o,m 
we introduce the notations 


(3.8) 

□ 

)„[X], 


■Kn{q) = {a I ai + ... + aN^ = q}, 
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and 


Nr, 

i=l 

Lemma 3.2. [X] admits the Wiener chaos expansion 




L2 


E E 

g—1 a£TT„iq) 



(3.9) 


where 

1 /* ^ 

a: jjRn \ ^ ^ / ■‘■■^ 

2=1 

„ TV 71 

^ y (JV -n) (^(2^^“)) n Har,+i{Ui)(p{Ui)du. 

2=1 

Proof. Consider 


(3.10) 

(3.11) 


F, 


(OjTTl)’ 


.[X] = 


0<,2i(A[;fV(i))det(Mx 


'(0,m) 


'(1) 


Aj'fiy 




{Afn-n} I ■'*■(2) 




ds. 


Take /i : M” —)• M and /2 : ” —)• M defined by 

fliv) = f 2 iu) = det (^Mj V{Nr,-n} (^( 2 ^^^) • (3-12) 

Then, /i e L‘^{R^,]Xl^^(i>{u)du), u e and /s G (j)iu)du), u G 

Write the Hermite expansions for /i ,/2 : 


where 


and 


71 


/l(«) 

II 

M 

■■■a„ HafiUi), 


>■" 

2=1 



Nn—n 

/2(ii) 

= E 

^an+i—aNn 


an+l,-" i^ATn 

2=1 


CL\ """OLyj 


Oi! 


-[ 

J Jw 


( 1 / 2 ) 


t2i(2 V(^) 




2=1 


1 

^T2+l! * * * • 


det (^Mi y{Nr,-n} (^(2(^^^) 

Nr, 

X Hai{Ui)(f{Ui)du. 

2 = 22+1 
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Then, with probability one, we have 
/l(y(l))/2(y(2)) = 
Re-arranging the sum gives 


Nn 


aTO+l-'-ajv, 




ANn 


i=l 


/i(r„|)A(yp,) = E E 

9=0 {ai}G7r„(g) 


Nr, 


‘a.n+l---a,Nr 




2 = 1 


= E E d^aHa{Ys)ds, 

q=0 a€TTn{q) 


yielding 


p CO 

eo,™)”W=7 e e 


To deduce the equality, write 


<3 


^Q= EE dZHa{Ys)ds. 

d(0,m)^ q=0SGnr,{q) 

Now note that the sequence {Aq}q^^ is Cauchy. To prove this, note first that 

- 2 

Mn+i—aN H^(Yg)ds 


WAqi-Aq^W"^ = E 


Q 2 


, , E E 


< mJ 


E 


' {o,my 


Q 2 


-\ 2 


E! E! daj^...a„da'n+l---aN„^a,iYs)ds 

q=Qi TTniq) 


Q 2 


< / E 

■^(0."^)" 5=Qi 


E! da^...aAo'n+l---aNnd^^As')ds 

Jn{q) 

where the last inequality here follows from the orthogonality of spaces TLq. 

Exploiting the independence of the compone nts of T, and applying a generalized 
Mehler ’s formula (see proof of pr oposition 2.1 in Estrade and Leon ( 2015l l and Lemma 
10.7 in Azai's and Wschebor ( 2009l lL we can bound the above expression by 


m 


Q 2 

“EE 

q=Qi T^n{q) 
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By convergence of the coefficients of the Hermite expansion, the above tends to zero when 
Qi,Q 2 increase, and so we have that is Cauchy. 

For the other faces, J G {J}\T° of dimension d = dim(J) < n, we have a slightly 
different expression for the coefficients dG in (13.9p . Recall (I2.42p . from which it follows, 
similarly to the above, that the corresponding integrands are given by 


h{v) = (vxj ) 1 




(3.13) 


hm = det (Mx, 
Thus, in terms of Y we have that 


(3.14) 


1 


JIR" 




„ N^—k 

X / _ det (^Mi^ R{jv„-n} H dia„+iiui)(p{ui)du, 

(3.15) 

where a! = oi! • • • and the second integral is evaluated over the Nk — k co¬ 
ordinates of u which appear in and the one coordinate that appears in 

^A-n-n) (aJT's). 

Lemma 3.3. F(o,m)"[A] admits the ITzener chaos expansion 

-^(O.m)" [Al] = '^q=i Ylnn(q) /(0,m)" dJ^{Yg)ds. 


□ 


(3.16) 


Alternatively, 


FiO,mr[X] = Y.W)- 

q=l 


(3.17) 


The variance, cr^, of F(^Q „i-^n[X] is given by 


— 


E E E 4-dga!6!R-(a,6), 

9=1 a&TTn{q) 6G7r„(q) 


(3.18) 


where G is given by 




a&iTn{q) 




symm ( 




20 




and the various coefficients are as follows: 


|det(A|J^^^)| (1/2) n 

d ,=—— n 

i=l 


i!(27r)’^/^ 


X 


„ Nn—n 

_^) det (^Mj JJ Ha^+i{ui)(j){ui)du, (3.19) 


ffi 


\a,b)=m'^f a'.bl ]J 




'(-m,m)" ^,^,>0 l<i,j<Nr, 

Yi dij=0,j 
j dij =bi 


l<k<n 




and 


^2 _ 
(Trr, — 


'Y 

q=l 


< = 


J] J] dsd-^albl R^{a,b), 

a^TTniq) b£TTniq) 


(3.21) 


Proof. By Lemma l3.II and Lemma 13.21 it suffices to establish the ffi convergence 

OO « CX) « 

Y Y Ha{Ys)ds^Y Y da H^{Ys)ds. 


q—0 a£xn(q) 


q—0 a£TT„{q) 


It is straightforward that hni —)• d^, where 


idet(A;;f))i-(V2) n 

d, = --n dda. (0) 

i=l 


d!(27r)”/^ 

X(iV -„) (^(^r^^)) (^(2^^“) n dIa„+i{Ui)(j){Ui)du. 


N„-n 


i=l 


We start by showing that (I3.16D is in ffi. By Fatou’s inequality 


E 


r r ~ 

Y Y da Hs{Yg)d. 

U=0ae7r„(g) d{0,mr 


1 2 


< lim E 

cr—>-0 


Q 


Y Y ds [ H^{Ys)ds 

g=0ag^4q) Y0,mr 


hmEV Y du [ Ha{Ys)ds 

-^0 J(0,m)- 


the last line following from orthogonality. 
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Adding some positive terms to the sum and then using Lemma FS-ll the above is bounded 


by 

2 

= E[F(^0,m)4X]f < oo. 

We introduce yet another shorthand notation, to be used for the remaining part of the 
current proof. 

/ Hs{Y^)ds and I, = V [ Hs{Yg)ds. 
ae^„(g) ae-n(g) 

With the above notation we have 

OO OO 

q=0 q=0 

Since lim dZ —)• dg, we have, for a fixed Q, 


- p 

—E / 

IT— >-0 _ n _ / N J(0.mV‘ 


HsiYs)ds 


lim 

(T—>-0 


Q 

E?? 

q=0 


L2 


Q _ Q _ 

E-^g = E ll^llra- 

q=0 q=0 


Moreover, since F(o,m)" [-^])-^(o ^ 'w® have 


(0,m)’ 


^[^]|| 


Now, 


E'.-E'J 

q=0 g=0 


< 


L2 


<9 




Q 


'WIL.' 


(3.22) 


E'»-E4" 

q=0 q=0 


+ 

oo 

E 4 

+ 

OO 

E 

L2 

Q — Q+1 

L2 

9— Q+1 


L2 


Given e > 0, we hrst choose Q' sufficiently large so that 

< s/3. 


E 4 

Q —Q^+1 


(3.23) 


L2 


Consequently, because of (|3.22l) . we can then choose a sufficiently small so that 


E 4 

9—Q^+1 


L2 


< e/3 and 


Q' 


Q' 


E 4 -E'? 

g=0 q=0 


L2 


< e/3. 


(3.24) 
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Since ||5^|| = 1, we relay on the fundamental relation for the Wiener chaos 

Nr, 


= 4 (symm(v9?“i 


(g)ajv„ 

^S,Nr, 


2 = 1 


(3.25) 


to write the expansion for [X], and then apply Fubini’s theorem for multiple Wiener 

integrals to arrive at 


« j-'in 

Q=0 7r„(<?) 

-eJ e 

9—0 \aGw„{g) ' ’ 


symm ( 


mNr, 

^S,Nr, 


ds 


9=1 


(3.26) 


with 


/™(Ai,--- ,A,)= ^ dJ 


deiTn{q) 


<{0,mY 


symm I 99^^ 




where Ai,... Ag G M, and M is as in (I2.19p . 

Next, we proceed to calculate the variance, Using (I3.16D and the orthogonality of 
T-Lq for different q, we have 




OO « 

V V dg / %{Ys)d^ 

9=0ae^49) 


H 2 




9-0 d£nr,(q) beTTn{q) 

By stationarity, this equals 

OO 

ZEE j 

9-0ae7r„(g)feg7r„(9) ^ ’ 

Then, a change in variables leads to 


E 


(0,m)^ 


Hn{Ys)HA% 


dsdu. 


(3.27) 


E 


Nr, 


Nr, 


n^a,(u,(o))n^6.m(^-^i) 


_ 2=1 


2 = 1 


dsdu. 


m 


OO « 

EE E J 

-n /.Nr* , , - 


E 


9-0 a£nr,(q) h£-nr,{q) 


{—m,mY 


Nr, 


Nr, 


ni?a,(ui(o))n^b,(^*(^)) 


_ 2=1 


2=1 


. l<k<n 
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When |a| = |6|, we have [see Proposition 2.2.1 in iNourdin and Peccatil (j2012l i] 


Nn 


Nr, 




i=l 


2=1 


dij > 0 l<i,j<Nn 


and zero otherwise. Thus, writing 

R^{a, b) = m" J 


n 1 


a' 6 ’ y TT 

f mm)" ■ ^ (daV- 

Y,i <^ij — 

Y. j dij = bi 


m 


the variance is given by 


dP, 


^2 _ 
— 


E E E d,d^,a\b\R^{a,b) = Y,<, 

9=1 aG7r„(g) 6e7r„(g) 9=1 


where u™ > 0 is given by 


< = 


J] J] dsd-.albl R^{a,b), 

Se7r„(g) b&iTr,{q) 


and we are done. 

Lemma 3.4. The coefficients in (j3.16p satisfy 

^ dlal<Cq^. 

d^-Kniq) 


Proof. The coefficients are 


|det(A<;f>)|-(‘/2) 




(27r)’^/2ai! • •-On! 




□ 


(3.28) 


Nn—n 

/ 2 (^ n ddar,+iisi)f)isi)ds, 

n 

2 = 1 

with /2 as defined in ()3.12p . It is straightforward to see that f 2 {^ is a polynomial of degree 
Nn — n + 1 and thus has a finite Hermite polynomial expansion. This means that all the 


1 


2n+l 


“iV„ 


an+1’ ■ ■ ■ ^Nn' J R^n — 
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terms da„+i-aN„ with any of the indexes a* > Nn — n + 1, z G {(n + 1),... , are zero. 
Setting 


c = |det(A;;;^OI 


( 1 / 2 )^|-( 1 / 2 ) 


X max 

^n + 1 


(4 




^72+1 • 




gives 


TTniq) ai,...,an£lTn{q) 


( HaM---HaA0) 

V ( 27 r)"'/ 2 ai! ■ ■ ■ aj- 


2 

ai! • • - On!. 


Using Imkeller et al. ( 1995l h Proposition 3, we have 


I ad 


< C and thus 


E 4+<C ^ 1<C,”, 

T^niq) ai,...,a„G7r„(ij) 

For other faces J G {J}\T°, by (I2.42p the contribution to the Euler integral of face J is 
given by the limits of expressions of the form 

J </’o-2i(VX|j(s)) det(V^X|j(s))X(s)l|i;vjf(^ 

which are bounded by 

(l + (det(v2X|j(s)))') (1 + (X|j(sl)2) ds. 

Although, the functions in the bound are slightly different to the corresponding functions 
in the previous development for the contribution of T°, the remainder of the argument is 
essentially the same, and so we shall not write out the details. □ 


3.2. Proof of Theorem \3.1\ 

As previously shown in Lemma f3.ll the sum corresponding to the interior critical points 
can be written as the limit of 

4omu[^]= / f>M^X{s))dei{V^X{s))X{s)ds. (3.29) 

By Lemmas 13.11 and 13.31 we have the Wiener chaos representation of (j3.29p 


FiO,m)AX]=YW)- 

q=l 


(3.30) 
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We, therefore, need to establish a CLT for 


m 


n/2 




m\ 

q ) 


q=l 


The CLT will follow immediately from Theorem 12.11 once we have checked that all the 
conditions of the theorem hold in our case. 

We start with Condition (c) oi Theorem EH Conditions (a) and (b) will be deduced 
from our previous results and the proof of (d). As far as (c) is concerned, note firstly that, 
by Lemma 13.31 we have 


/-(Ai,... ,A,)= ^ [ 

a^iq) 


symm I 


(g)ajv^ 

^S,Nr. 


ds. 


where 


symm ® ® ) = ^ Z] ® ® ’''' ’ ^-(-?)) ' 


Since |a| = q, the inner sum can be written as 

Nr, 






(3.31) 


with the appropriate coefficients i 9 =i’ = 0, whenever YliJi / 

q. Take C{q) = Using the above notation, we write 


aG7Vn{q) b^TTniq) 

Z! {}Ps,jx <8> • • • <8> q^sjq) ®r {^u,ki 


'(0,m)2 


' ^u,kq) dsdu. 


31,--^Jq 

ki ,...,,kq 


The following is true (see Section 12.5p 

® ® ^xjq^ 'Sir (^'Py,ki S ' ' ' S 'Py,kq'} 


'Py,ki)sj 


1=1 


fx,jr + l S''' 'Pxjjq S fy^kr+1 S ' ' ' ^y,kq 


d^jiM {x y)'px,jr+l S ■ ■ ■ 'Px,jq S l^y^kr+l S ' ' ' 'Py,kq 1 


1 = 1 
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and 


1=1 


(v^xjr+l ® ■ ■ ■ ‘^x,jq ® ‘^y,kr + l ® ' ' ' '^y,kqj 9^-W,V+l ® ' ' ' ‘^W,iq ® ¥^2,mr+l ® ■ ■ ■ '^Z,mq')^^2q 

q—r q—r 

= '[l^jl,ilix-w)'[lKki,miiy- ^■ 

1=1 

We have 

r 

n {x-y)< ip^{x - y), 

and 

q—r 

Y[^ji,iiix - w)Y[Kki,mi{y- z) < - z). 

1=1 1=1 

Since the number of summands in (I3.3ip is less then [Nn^, we have that 


1=1 


q—r 


I E 4] iiNnfcw)fz{t), 

,aG7r(g) 


where 


Z{m) = f {x — y)'ilj'^{w — z)ili'^ ^{x — w)ijj'^ ^{y — z)dxdydwdz. 

J 


Next, using ip^{x — y)iji^~^{x — w) < 'ip‘^{x — y) + '4’‘^{x — w) we have 

Z{m) < Zi{m) + Z 2 {m), 

where 

Zi{m) = / {w — z)il}'^{x — {z — y)dxdydwdz, 

Z 2 {m) = j {w — z)il)^{x — y)il}^~'^ {z — y)dxdydwdz. 

J(0,m)4" 

Integrating with respect to x and using 


/ 'iIj^{x — w)dx < / ijj‘^{v)d' 


'v < OO, 
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then integrating with respect to the remaining coordinates and using 
f {w — {z — y) * dydwdz = 

— y)dwdy < f ip^ *{v)di 


we get 




Y, 4 ) {{NnycfZ{m)<oo, 

, ae7r(q) 


from which it follows that 


Consequently 




< 


c 


m" 


^ ^ ^ ^ Yn^oo 

Thus, we have established that condition (c) of Theorem IQ is satisfied. 
We now turn to the condition (d). We have to show that 


supm " Y 


m\\2 


m>l 


Q —Q+1 




0 . 


The expression i® the variance of the tail of the Wiener chaos 

expansion of We need to show that it converges uniformly to zero. We 

have already developed this expansion in the previous results, with the only difference that 
now we have a normalization of m~^, 


1~Q aS7r„(g) b£iT„{q) 


(3.32) 


X / E 

J 


'N„ 


Nn 




. 2=1 


2=1 


n (1 - '^)du. 


l<k<n 


To show the uniform convergence we write an upper bound C{Q) which is ind ependent 
of m and vanishes as Q —)• oo. To construct such a bound we use Lemma 1 in ArconesI 
(jl994l L which, for completeness, we reproduce. 
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Lemma 3.5 (jArconesI ([199J)). Let V and W be two zero-mean Gaussian random vectors 

^ M have Hermite rank r. (i.e. 


on 


and assume that "EViVj = 'KWiWj = 5ij. Let h 


the lowest degree polynomial appearing in its Hermite expansion has degree r.) Write ip* 
for the supremum of the sum of the rows or columns of the covariance matrix Cov{V, W), 
and assume that i/’* < 1. Then 

E[h{V) - Eh{V)][h{W) - Eh{W)] < iflEh^iV). 

Returning to the proof of Theorem 13.11 we now apply this lemma with 


R = (Fi(0),... , TwJO)), W = (yi(z7),..., YjvJi7)), 


and h : —>■ IR given by 

Nr, 

Hx) = '[lHa,{Xi). 

i=l 

It is easy to check that 


V’* < K‘^iP‘^{t) and Eh^ = ^ d| a!. 

T^nlq) 

Furthermore, since <? > 0, IEO^l = 0, and for |t| large enough, by the assumption 

on ip{f), we have that K^X{t) < 1. 

We now choose arbitrary s € M'*' and split the integral over two domains 


E 




Nr, 


Nr, 




_ 2=1 


i=l 


l<k<n 


n 


E 




Nr, 


Nr, 




i=l 


2=1 


n b 


l<k<n 


\Jfk\ 

m 


did 


+ / E 

J{—m,mY\Rg (s) 


Nr,. 


Nr,. 




. 2=1 


2=1 


n (' 


\^k\ 

m 


dv. 

(3.33) 


Here Hq(s} is n-dimensional cube of side length s, centered at the origin. 
For the first term corresponding to the integral over i?g(s), we write 



Nr, Nr, 

dr,\{Hr,^{Ym) E d^HHbAUv)) 

a&-Kr,{q) *=1 beirr,(q) 


n 

l<k<n 



did. 
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Lemma 13.31 implies that the above sum is finite for all Q, and so it converges to zero 
as Q —)• oo. Regarding the second term, and reintroducing the summation from (13.321) . 
consider 


E 


E 


q=QA-m,m)^\R^{s) 


Nr, 




_ iVn 

b^TTniq) 


X n L 

l<fc<n ^ 



di7. (3.34) 


By Lemma 13.51 we have 


E 


Nr, Nr, 

^ d.HHa^O)) 

a&-Kr,{q) *=1 b&irr,{q) 


< Y, dld\, 

■^niq) 


so we can bound the second integral by 


E / 

g^QJi-m,mY\K^{s) 




l<A:<n 




< 


oo 

sE 


K'^'ip^{V) Y] d^aldu. 


By Lemma [R4l we have d|a! < Cq^, so that 

a&TTr,{q) 

[ I 4 a! I di7< [ K<i'ilj<i{i7)Cq^diJ. 

J{-m,mY\R^{s) \'Kr,iq) j J i-m,mY\R^ (s) 

Due to the assumption that V’(^) 0, we have, for s € M'*' large enough, that 

1 




R’ 


(3.35) 


which leads to the bound 


Y Cq^ - / mdi 


q —Q+1 


Since Ke < 1 and 


/ ii{R)d\ 

JMP 


V < OO, 
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the integral in (j3.34l) converges to zero uniformly in m. Overall, we conclude that m —)• 

Q—>-oo 

0 uniformly in m, which establishes that Condition (d) is satisfied. 

To show (a) and (b), we note hat we have already demonstrated that for sufficiently 
large s the integral over (—m, (s) converges to zero. Then, 


R^{d,b) ^ R{a,h) = / a\b\ 


n 


E 

dij>0 

dij=aj 

dij=bi 








leading to 


OO 

2 2 ^ 

(7rfi ^ 
m^oo 

q=l 




where 


(3.36) 


Uq= ^ ds d^a\b\R{d,b). 

0,G7Tn{q) bGTTn{q) 

This completes the proof of Theorem 13.11 


□ 


4. The mean value of upper Euler integral 


In this section, we discuss the mean value of the Euler integral. As shown bv iBobrowski and Borman 


(120121 ) , the mean value of the Euler integral of Gaussian random field scales not by the vol¬ 
ume of the domain of integration, as one would expect, but according to a one-dimensional 
measure of the domain. Specifically, 


E[Tm[X]] = - 


^f(M) 


(4.1) 


where C^{M) is the first Lipschitz-Killing cur vature of M, as e v aluat ed with respect to 


the metric induced by the random field X. (cf. Adler and Tavlor ( 200?! ) for definitions.) 

In this section, we re-establish this result by direct evaluation of the mean value through 
the Wiener chaos decomposition of the Euler integral. To do so, we make use of the next 
proposition, which can be proven using symmetry considerations. 


Proposition 4.1. Let X be tame Gaussian on M”, and n > 1. Then, 

E [det(V2A)A:] = 0. 


(4.2) 
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Recall that in proving the CLT of the previous section we concentrated only on critical 
points in the interior of the parameter space which contributed to the Euler integral. Now, 
however, we need to consider all such points, since we are looking at an un-normalised 
mean, rather than an asymptotic limit. 

What is now interesting, and very different to what we saw before, is that the chaos 
approach shows that none of the faces of dimension different from one can contribute to 
including the interior face. This gives, from this angle at least, some new 
intuition into the Bobrowski-Borman result. 

To justify this claim, note that in our chaos expansions the mean values of random 
variables in a chaos of order greater than zero vanish, so that possible contribution to 
mean values may come only from the zeroth chaos. This, however, is characterized by the 
coefficients do-'-olJ- Let {d^}s be the coefficients of the chaos of some face J. In general, 

each in the Wiener chaos expansion for J of dimension k factorizes as - , 

where 

n 

y.WHai{yi)(j){vi)dv, (4.3) 

i=l 

Nk-k 

X Ha„+iiui)(j){ui)du. (4.4) 

i=l 


d^^^ I T = lim — 
cr-j-o a! 


Vij 


d 


( 2 ) 


L 




det (Mj 


Note, we have that ~ IL[det(V^X|j)X|j]. Thus, by Proposition 14.11 the contri¬ 

bution to the mean value of the faces J, dim J > 1 is zero. Since the underlying field is 
centered it is obvious that the zero dimensional faces, the vertexes, do not contribute to 
the mean as well. Overall, the conclusion is that only the edges contribute to the mean 
value of Euler integral. 

To see what this implies in a simple example, take M = with additional assumption 
of isotropy. Then 


p) 

^o,...,o 


= E[ 


d‘^X{s) 

dsiSi 


xm = - 


(9^/9(0) 

Ondn 


— -A 2 , 


(4.5) 


where A 2 is the second spectral moment of X. 
To calculate g, we consider 


Pi{VX\j{s)))det{V^X\j{s))X{s)p^^xi^^^j)m- ( 4 - 6 ) 


32 




Having summed over all parallel edges and taking expectations, we can eliminate the 
indicator term l{(vx(s),r;j)>o}) since by stationary we translate everything to the same 
range over (0, m) which yields 

= 1- (4-7) 

{Ji\Ji parallel to J} 


Consequently, by (I3.19p . 


,(i) ^ 

(27r)V2 ’ 


(4.8) 


and since the set of the edges of can be split into n families of parallel edges, we finally 
have 


= 4 '^ X 4 ^ = 


1 ^ 2 ] X m 

(27r)V2 


(4.9) 


This is precisely (14.ip for this case. 
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